Abstract. We analyze new high-resolution Chandra X- 
ray images of the Wolf-Rayet binary system WR147. This 
system contains a WN8 star with an early-type compan- 
ion located 0.6" to its north, and is the only known early- 
type binary with a separation on the sky large enough 
for the wind-wind collision between the stars to currently 
be resolved at X-ray energies. The 5 ksec Chandra HRC-I 
image provides the first direct evidence for spatially ex- 
tended X-ray emission in an early-type binary system. 
- ■ ■ The X-ray emission peaks close to the position of the 
. radio bow shock and north of the WN8 star. A deeper 
X-ray image is needed to accurately determine the de- 
C\| . gree of spatial extension, to exactly align the X-ray and 
$_i ' optical/radio frames, and to determine whether part of 
Oh. the detected X-ray emission arises in the individual stellar 
winds. Simulated X-ray images of the wind-wind collision 
have a FWHM consistent with the data, and maximum 
likelihood fits suggest that a deeper observation may also 
constrain the inclination and wind momentum ratio of this 
system. However, as the WR wind dominates the collid- 
ing wind X-ray emission it appears unlikely that M qb and 
^^oooB can be separately determined from X-ray observa- 
tions. We also note an inconsistency between numerical 
and analytical estimates of the X-ray luminosity ratio of 
,: the stronger and weaker wind components, and conclude 
' that the analytical results are in error. 
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1. Introduction 

Massive, early-type binaries with powerful stellar winds 
can generate a complex region of shock-heated plasma 
with temperatures in excess of 10^ K when these winds 
coUide (Prilutskii & Usov |l976| ; Cherepashchuk |l976D . X- 
ray observations provide a direct probe of the conditions 
within the collision zone and also of the unshocked atten- 
uating material along the line of sight in the system (e.g. 
Steve ns et al. |l992| ; Pollock et al. |l999t Pittard & Stevens 
1997). Single early- type and Wolf-Rayet (WR) stars are 



also known to be X-ray emitters. The most viable model 
currently proposed is that the 10^ — 10^ K X-rays are 
generated by shocks in the wind resulting from the unsta- 
ble nature of the line-driven acceleration (e.g. Owocki et 
al. 1988; Feldmeier et al. 1997). Direct spectroscopic evi- 



dence of clu mps/sh ocks in WR winds was first noted by 
Moffat et al. ( |l988D . 

Emission from colliding winds is often observable as 
a strong high-temperature X-ray excess above that ex- 
pected from the individual stars, together with phase- 
locked orbital variability resulting from the changing line 
of sight into the system (e.g. Pollock 1987; Chlebowski 
1989| ; Williams et al. |l99(]| ; Chlebowski & Garmany |1991 ; 
Corcoran 1996). Recent Chandra observations have re- 
vealed spectral features which are noticeably different to 
those from single early-type stars, including strong forbid- 



den line emission (Corcoran et al. 2001 ; Pollock et al. 2002 
Skinner et al. pOOll ). In all previous X-ray observations, 
the colliding wind emission has remained spatially unre- 
solved, and contaminated by the emission from the indi- 
vidual winds. 
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Non-thermal radio emission from WR stars is also a 
good indicator of wind-wind interaction (e.g. Dougherty 
& Williams 200C). Such emission is thought to be syn- 
chrotron radiation arising from electrons accelerated to 
relativistic velocities in the wind-wind collision zone (Eich- 
ler & Usov 1993 ). The presence of nonthermal radio 
emission fr om W R 147 (AS 431) was re ported by Cail- 
lault et al. (1985) and Abbott et al. (1986). Later observa- 
tions with MERLIN at 5 GHz resolved two radio sources 
( WR 14 7S and WR 147N) separated by ~ 0.6"(Moran et 
al. 1989). The southern thermal source (WR 147S) was 
coincident with the WR star, and the northern source 
(WR 147N) was proposed to be responsible for the ob- 
served nonthermal emission. Independent observations 
with the VLA confirmed this hypothesis (Churchwell et 
al. |1992| ; Skinner et al. |1999D . Moran et al. ( |1989D further 
suggested that WR 147 could have an unseen companion 
associated with the northern radio source. 

WR 147 attracted the attention of many astronomers 
when radio, infrared and optical observations at high spa- 
tial resolution provided unambiguous evidence of its bi- 
narity (WiUiams et al. 1997, hereafter W97; Niemela et 
al. |1998| , hereafter N98). New, higher resolution MERLIN 
observations by W97 revealed that at 5 GHz both sources 
are elongated with sizes of ~ 170 x 253 mas (WR 147S) 
and 267 x 79 mas (WR 147N), and that they had a sep- 
aration of 575 ± 15 mas. The UKIRT infrared and HST 
optical images revealed for the first time a faint companion 
star located to the north of the nonthermal source associ- 
ated with WR 147N. Crucially, the relative positions of the 
sources revealed that the companion star, although close 
to the nonthermal radio emission, was located slightly (« 
60 mas) more distant from the WR star. The resulting 
picture is thus consistent with the nonthermal emission 
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arising in a colliding winds shock between the two stars, 
with the wind of the WN8 star dominating. 

The spectral type of the companion is currently some- 
what uncertain, as the excessive reddening towards this 
object {Ay ~ 12 mag; Churchwell et al. 1992) makes it 



2. HRC-I Analysis 



hard to get a (blue) classification spectrum. Based on ram 
pressure arguments, a companion with a late-0 or early- 
B spectral type is expected to have a sufficiently strong 
stellar wind to produce the bow-shock at its observed po- 
sition|^. Although the mass-loss rate and wind velocity are 
not known for the OB star, the wind momentum ratio can 
be directly determined from the geometry, specifically by 
comparison of the position of the nonthermal source with 
those of the stars (i.e. the ratio ro^/D, cf. Usov 1992| ). 
Estimates for the wind momentum ratio, 

(MWoo)oB 



V 



{MVoo)wN 

range from 0.011 



+0.016 
-0.009 



(W97) to 0.028 



^-Vo?7 (N98). 



(1) 
This 



value can then be used to estimate the wind momentum 
of the OB star if we know that of the WN star (see Sec- 



tion 3.1 ). The inclination angle can also be estimated from 
the observed geometr y: W 97 deduced i — 41°, and Con- 
treras & Rodriguez ( |1999| , hereafter CR99) determined 
z = 45 ± 15° from VLA observations. 

With WR 147 revealed as a colliding winds source it 
was expected that most of the X-ray emission could be 
from the wind-wind collision. A recent observation with 



ASCA (Skinner et al. 1999) determined the emission to be 
thermal (although multi-temperature) with the dominant 
component at kT « 1 keV. This is slightly harder than 
the characteristic temperature from single early-type stars 
(« 0.5 keV). A high absorption column (TVh ^ 10^^ cm~^) 
in agreement with estimates based on the visual extinction 
was also confirmed. The intrinsic luminosity in the 0.5 — 
10 keV band was estimated as Lx = 3.5 x lO'^^ ergs~^ , 
which gives log Lx/Lboi ~ —6.5, 0.5 dex greater than the 
canonical value of —7 from single early- type stars. 

WR 147 is the second closest WR star (D « 650 pc; 
Morris et al. 2000) and one of a small number of WR sys- 
tems where the companion stars have been spatially re- 
solved (other examples being WR 86 and WR 146 - N98). 
Of these, only WR 147 has a sufficiently great separa- 
tion between the stars for us to attempt to spatially re- 
solve the wind-wind interaction region at X-ray energies. 
This in turn is currently only possible with the recently 
launched Chandra X-ray Observatory, which combines a 
superb mirror with sub-arcsecond resolution (van Spey- 



broeck et al. 1997) and accurate, stable pointing. In this 
paper we report on the analysis of a 5 ksec observation of 
WR 147 with the High Resolution Camera (HRC-I). 



^ The latest attempt (Lepine et al. |200l[ ) based on STIS/HST 
resolved spectra yields an 05- 71(f) type, which seems a little 
luminous for the observed magnitude difference. 



Chandra (Weisskopf et al. |1996| ) obs erved WR 147 with 
the HRC-I instrument (Murray et al. |l997| ) on 2000 July 
10 12:11 UT to 14:09 UT for a total exposure time of 
4.87 ksec. The HRC-I is a microchannel plate imager with 
excellent spatial and temporal resolution but poor energy 
resolution. Each photon detected by the HRC-I is time 
and position tagged. The size of the electronic readout 
is 6.429 fim (0.13175", referred to as 'one pixel') which 
oversamples the point-spread- function (PSF). All Chan- 
dra data which have undergone reprocessing have bene- 
fited from improvements made to the crossed-grid charge 
detector signals. This has reduced the 50 per cent encir- 
cled energ y radi us from 4.0 to 3.2 pixels (0.527" to 0.422"; 
Juda et al. pOOO ). Softer photons are brought to a slightly 
tighter focus than harder energies. Our analysis is based 
on the pipeline processed data from the Revl Chandra 
X-ray Center (CXC) reprocessing. 

The CXC has measured the on-orbit performance of 
the Pointing Control and Aspect Determination (PCAD) 
system on Chandra (see Sec. 5.4 and Table 5.1 of the Pro- 
poser's Observatory Guide, hereafter POG), and deter- 
mined that the standard processing is capable of placing a 
reconstructed X-ray image on the celestial sphere to an ac- 
curacy of 0.76"(RMS) radius. The image was dithered dur- 
ing the observation and these effects were removed during 
post-processing. We also checked that pointing transients 
had damped out by the time that the actual observation 
began, and that the pointing remained stable during the 
observation. 

A total of 115420 events were noted after the standard 
level 2 processing, which yields a count rate across the 
entire detector of 23.7 cts s~^. This implies a background 
rate of slightly less than 10~^ cts s^^ arcsec"^. We there- 
fore expect only 1 background count within the source 
area of WR 147 shown in Fig. |l|. Running the source de- 
tection software yielded one reliable source in addition to 
WR 147, "Source 2" in Table |. We did not find an optical 
counterpart to this source in the Tycho- 2 catalogue from 
the Hipparcos mission (Hog et al. 200C), but a counter- 
part was found on images from the Palomar Sky Surveys 
(POSS). Its position has been measured on red and blue 
plates of the (original) POSS-I survey with the PMMma- 
chine for the USNO-A2.0 catalogue (Monet et al. |l999D 
to be 20 36 30.100, -f40 21 23.32 (J2000). The offsets 
in RA and Dec from the Chandra detection are -fOs.018 
and -fO.22", close to the estimated astrometric precision 
(0.3") of the USNO-A2.0 catalogue. The I magnitude of 
the coincident source is ~ 12, which is slightly too faint 
for Tycho. 

We also measured the star's coordinates from Super- 



Cosmos (e.g. Hambly et al. 1998) scans of glass copies of 



the POSS-II survey. It appears on the overlap region of two 
red plates giving coordinates 20 36 30.058, -t-40 21 22.97 
and 20 36 30.074, +40 21 22.91 (J2000). These posi- 
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Fig. 1. Raw X-ray image of WR 147 containing 148 counts. There are 3 pixels towards the centre of the image which 
have the maximum of 4 counts per pixel. Only 1 background count is expected in this image - all of the others are 
from WR 147. 



tions are independent of each other, coming from reduc- 
tions using different sets of Tycho2 stars centred about 
5 degrees apart. The mean offsets from these coordinates 
to the Chandra Source 2 are -0.014s and -0.16", in the 
opposite sense to the POSS-I offsets. The sets of stellar 
coordinates from the POSS-I and POSS-II plates arc con- 
sistent within the uncertainties of their measurement but 
the difference (0.54") might also reflect proper motion of 
the star between the epochs of the POSS-I and POSS-II 
survey observations. 



Finally a search in the 2MASS catalogue also revealed 
a possible match for the second source, although the po- 
sitional accuracy is too poor to help to register the HRC- 
I image. In any event, because we have only one X-ray 
source additional to WR 147 in the field, which could be 
a linc-of-sight coincidence with the star measured off the 
POSS, we will not use these offsets to refine the align- 
ment of the X-ray and optical/radio positions, and must 
accept some astrometric uncertainty in our analysis of the 
Chandra data. 
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Table 1. List of detected sources from the wavdetect algorithm. The extended emission at the position of WR 147 
is detected with the greatest significance, and the quoted positions are from an unblocked image. All uncertainties, 
including the count rates, are determined from the wavelet analysis. 



Source 


RA (J2000) 


Dec (J2000) 


Net counts 


WR 147 


20 36 43.637 ± 0.003s 


+40 21 08.10 ±0.04 


148 + 12 


Source 2 


20 36 30.082 ± 0.010s 


+40 21 23.05 + 0.10 


19 + 6 



The accuracy of the radio positions of WR 147 and 
their tie to the optical frame comes via the phase refer- 
ence source, and the uncertainty is not much more than 
10 mas in each coordinate. On the other hand, the uncer- 
tainty in the X-ray positions relative to the optical frame 
is a combination of the source centroid on the detector 
(as listed in Table |^), and of the tie of the X-ray to the 
optical frame. It is this last factor which dominates the 
uncertainty in the relative radio/X-ray positioning. As al- 
ready noted, the standard processing is accurate in a sta- 
tistical sense to 0.76"(RMS). However, the Chandra and 
SuperCosmos-POSS-lI positions of source 2 suggests that 
the positioning of the X-ray frame is accurate to -0.014s 
and -0.16", which is consistent with the uncertainties in 
the position of the source centroids on the X-ray detector 
and also comparable to the SuperCosmos uncertainties. If 
we were to incorporate this shift, the X-ray peak would 
lie on the western edge of the radio bow shock. Hence the 
intrinsic astrometry appears to be much more accurate for 
our observation than indicated by the RMS uncertainty, 
and we can confidently state that the X-ray peak is defi- 
nitely not cospatial with the WR star. 

However, we would need to be able to align the X-ray 
and optical frames to a positional accuracy of better than 
0.05"(< 0.4 pixels) to distinguish between the position of 
the radio bow shock or the companion star. This is clearly 
beyond the capability of our current dataset. Therefore in 
Fig. I we have added an offset (-0.00965s, -0.1387") to 
the actual X-ray data so that the peak emission matches 
up with the centre of the nonthermal radio emission from 
the bowshock. This offset is well within the nominal ac- 
curacy of the Chandra aspect solution. In summary the 
X-ray peak is certainly not cospatial with the WR star, 
but could be cospatial with either the nonthermal radio 
peak or the northern optical component. 

Another crucial question concerns whether the de- 
tected emission is consistent with an extended object. 
To address this point we constructed several independent 
tests. We first extracted a PSF at the position of WR 147 
using monochromatic source photons of energy 1.5 keV. 
This is a sensible approximation to the actual spectrum 
which is sharply peaked at « 1.5 keV (Skinner et al. 19991 ). 
We then extracted radially averaged surface brightness 
plots of the PSF and WR 147, and found that the full 
width half maximum (FWHM) of the former was ~ 0.45" 
(in good agreement with the 0.4" quoted in the POG), 
compared to « 1.0" for the WR 147 data. 



In a second test we fitted a Gaussian function to the 
PSF and the HRG-I data of WR 147, and used a maxi- 
mum likelihood method to determine the optimal FWHM. 
Despite experiencing complications in our fit to the PSF 
data due to its extended wings, we again found that the 
FWHM of the WR 147 data (« 0.85") was significantly 
larger than that of the PSF (w 0.5"). Finally, in a third 
test we investigated how the likelihood varied if we as- 
sumed that the source was a circular region of constant 
surface brightness and adjusted its radius. When fitting 
to the PSF, we found that the likelihood declined rapidly 
as the radius was increased. In contrast, the likelihood was 
almost constant for a source radius between and 5 pix- 
els when fitting to the WR 147 data. (i.e. for a circular 
source up to a radius of 0.66"). This contrast in behaviour 
again supports the inference that WR 147 is extended. We 
therefore conclude that WR 147 is resolved at X-ray en- 
ergies, and note that to the best of our knowledge this is 
the only stellar system where this is the case. 

While the HRG-I has poo r spe ctral resolution, previ- 
ous ASCA data (Skinner et al. 1999| ) determined that there 
were essentially no X-ray photons from WR 147 with en- 
ergies between 0.4-1.0 keV. Since the extinction by dust 
towards WR 147 is very high, we conclude that there es- 
sentially all X-ray photons with energies below 1.0 keV 
are absorbed. It is therefore highly likely that (at least 
the majority of) the emission detected in the HRC-I is 
from the wind- wind collision zone, and we again note con- 
sistency with our conclusion that we are seeing extended 
emission from WR 147. 



3. Theoretical X-ray Emission Maps 

Based on the above conclusions, we now explore the fea- 
sibility of using X-ray imaging to constrain some of the 
fundamental parameters of this system. This involves com- 
puting a grid of synthetic emission maps from hydro- 
dynamical calculations of the wind-wind collision. Since 
there is some small possibility that each star is also con- 
tributing to the observed emission, we also attempt to 
determine the relative luminosities of each of the three 
possible components (i.e. each star plus the wind-wind 
collision) . 
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Fig. 3. Hydrodynamic simulation of the colliding winds 
in WR 147. The WN8 star is located at the bottom left 
of the plot at position (0,0), with the companion at po- 
sition (8.8 X 10^^,0). Shown is a gray-scale of the den- 
sity (logiop g cm~^). It is clear that the collision region 
is globally stable. The parameters used in this model 
(cw_2) are: M wn = 1-5 x 10"^ Mgyr"!, M qb = 
2.8 X 10""^ Moyr^i, Woown = 950 kms~S Wooob = 
1000 kms-1 . 

3.1. Estimates of the Wind Parameters 

To estimate the X-ray luminosity from each component, 
we first need to obtain values for some essential parame- 
ters, such as the mass-loss rates of the stars. 

On account of its proximity, WR 147 has been ex- 
tensively studied over different energy bands, and the 
wind parameters of the WN8 star are well determined. 



Tight constraints were derived by Morris et al. (2000) us- 
ing ground-based optical and near-infrared observations, 
combined with high resolution space-based ISO obser- 
vations. They obtained a dumping-corrected mass-loss 
rate of Mwn = 1-5 x lO~^M0yr~^ (with filling fac- 
tor /(radio) ~ 10%), substantially lower than derived 
from optical or radio observations under assumptions of 
homogeneity. The terminal wind velocity, obtained from 
helium recombination lines and fine structure transitions 
in the infrared, was also revised downwards to Woown 
950 kms~^ . The abundances of the WN8 wind (which af- 



fect our calculation of the X-ray emission from the wind- 
wind collision in Secti on |3.3| ) are also fairly well deter- 
mined (cf. Morris et al. 2000| ). From information contained 
in their paper we have adopted the following relative abun- 
dance by mass (sum = 1): H = 0.09, He = 0.9, C = 10"^, 
N = 0.01, O = 8x10-*, Nc = 1.2xl0"^ Mg = 4.8xl0"^ 
Si = 10-3, S = 2.1 X 10""*, Ar = 2.4 x 10-^ Ca = 
4.7 X 10-^ Fe = 2 x 10-^ Ni = 8.8 x 10"*^. For the OB 
companion we assume solar abundances. 

As for almost every colliding winds system, the wind 
parameters of the companion are not well determined. It 
is not even possible to place a reliable upper limit on the 
mass-loss rate of the wind of the companion star from the 
non-detection of flux in the high resolution radio images, 
since the position of the star is on the edge of the nonther- 
mal source where the flux gradient is very steep. Further- 
more, the values of rj derive d by W97 and N98 depend on 
the value of tob (cf. Usov 1992), which is the difference 
of two rather similar separations, each with observational 
errors, and hence very uncertain. We therefore choose to 
average the results of W97 and N98 in the observational 
frame, and determine a projected stellar separation of 
639 mas. Then using the projected radio separation de- 
duced by MERLIN (575 mas), we obtain = 0.0121°;°°^, 
independent of inclination. 

Given this wide range, we investigate the X-ray emis- 
sion properties from a series of colliding winds models with 
varying r] and inclination, i, as listed in Table ^j. To reduce 
the number of free parameters we fix the terminal veloc- 
ity of the companion star at 1000 kms"^ , although there 
is plenty of scope for a higher velocity. With all of the 
other parameters fixed, this leads to the mass-loss rate of 
the OB star varying directly with the chosen value of rj. 
As it is possible that the X-ray emission from the wind- 
wind collision is somewhat contaminated by emission from 
the individual winds, model X-ray luminosities lower than 
observed are to be preferred, and we have weighted our 
parameter space accordingly. 

3.2. X-ray Emission from the Stars 

We first attempt to estimate the X-ray luminosity from the 
WN and OB stars. There is a broad consensus that X-ray 
emission from early-type stars is generated from shocks 
withi n the stellar wind (cf. Owocki et al. 1988; Feldmeier 
et al. 1997). However, it is very difficult to make reliable 
estimates of the intrinsic emission from each star as the 
variation in the observed luminosity between similar ob- 
jects can be very large. Despite this, the following calcu- 
lations give a rough impression of what we might expect 
for the intrinsic luminosity from each star. For OB stars, 
Sciortino et al. (1990) note a strong correlation between 
the intrinsic X-ray emission and the wind momentum rate: 



log = 0.53l[;;f^ log(Mfoo) + 17.8^ 



+3.5 
2.7 



(2) 



Substitution of the values of M and Voo from Table g 
gives Lx ~ 10'^^ erg s"^ (although with an uncertainty of 
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±3 orders of magnitude) . Taking the central values on the 
constants in Eq. || (and assuming a 0.01-10.0 keV band) 
we find that this is between 5 and 35% of the intrinsic 
colliding winds emission from our models (see Section pT3| ) . 

Estimates of the intrinsic X-ray luminosity from WN8 
stars have been made by Wessolowski ( |1996D and Pol- 
lock ( |19871 ), and again the variance between objects can 
be large. This is partially caused by large uncertainties 
due to the small exposure durations that many of these 
stars were observed for. Therefore, perhaps the most ac- 
curate estimate can be made from the i? 05^47" observa- 
tions of WR 16 (HD 86161; WN8h, v=8.44) and WR 40 
(HD 96548; WN8h, v=7.85), neither of which were con- 
vincingly detected during WR 16's 7554s PSPC exposure 
or WR 40 's 45386s HRI exposure, giving < lO^iergs^i 
for WN8 stars. Thus the OB companion is more likely 
than the WR star to contaminate the emission from the 
wind-wind collision zone. 



3.3. X-ray Emission from the Wind- Wind Collision 

In this section we determine some of the basic properties 
of the wind collision region, using both analytical and nu- 
merical methods. We use the results from previous obser- 
vations as a starting point, and therefore assume that the 
wind momentum ratio lies between 0.005 — 0.02 and the 
orbital inclination between 30-60° (cf. W97, N98, CR99). 
The full range of computed models is listed in Table |[ 

First we use an analytical approximation (cf. Eichlcr 
& Usov 1993) to estimate the half-opening angle of the 
contact discontinuity, measured from the line between the 
secondary star and the shock apex. For the models in Ta- 
ble I we find that ~ 31° (models cw_l - cw_3), 6 ~ 26° 
(models cw_4 — cw_6), and ~ 19° (models cw_7 — cw_9). 

The combined kinetic power of the two winds for 
the parameters in Table is approximately constant at 



To obtain accurate estimates of the X-ray emission 
from the wind-wind collision, and to generate synthetic 
emission maps with which to compare to the Chandra 
data (see Section ^^), we have computed hydrodynam- 
ical models of the wind collision, calculated using VH-1 



4.3 X 10 ergs" (as it is dominated by the WN8 wind 
for which we assume fixed parameters). An indication of 
which wind dominates the X-ray emission can be found 
(to first order) by evaluating the characteristic cooling pa- 
rameter, Xi of each wind. To do this we simply make use 
of Eq. 8 in Stevens et al. ( |1992| ), while noting that WN8 
abundances lead to a similar overall emissivity as solar 
abundances. Values for xwn and xob are listed in Table ||, 
and since xwn and xob are always ^ 1, the emission from 
the shocked winds is close to the adiabatic limit. In such 
circumstances, the denser wind dominates the emission 
(cf. Myasnikov & Zhekov |l99|), and we indeed find from 
hydrodynamical models that icw.WN ~ 20Lcw,ob- In this 
regard, however, the results from hydrodynamical models 
are in conflict with estimates from the relevant equations 
in Usov ( 1992| ), which predict that Lcw,wn ~ icw,OB- This 
is a worrying revelation as these equations are widely used 
in the literature, and we investig ate th is discrepancy in a 
further paper (Pittard & Stevens 2002). 



(Blondin et al. 1990), a Lagrangian-remap version of the 
third-order accurate Piecewise Parabolic Method (PPM; 
Colella & Woodward 1984). The stellar winds were mod- 
elled as ideal gases with adiabatic index 7 = 5/3. Although 
the wind collision is largely adiabatic, for completeness 
we included radiative cooling via the method of opera- 
tor splitting. The cooling curve for the temperature range 
4.0 < log T < 9.0 was generated usi ng th e Raymond- 
Smith plasma code (Raymond & Smith 1977), which spec- 
ifies an effectively thin plasma in ionization equilibrium. 

The simulations were calculated with the assumption 
of cylindrical symmetry - the orbital period of WR 147 
is likely to be of order thousands of years, so all orbital 
effects are negligible. The large orbital separation also im- 
plies that the winds are likely to collide at very near their 
terminal velocities. We therefore assume that we can treat 
the winds as being instantaneously accelerated to their 
terminal velocities at the stellar surface of each star, and 
do not consider radiative driving effects. For simplicity in 
this first investigation, we also assume that the post-shock 
ion and electron temperatures equalize very rapidly. While 
the equilibration timescale can be significant compared to 
the flow timescale for wide systems (e.g. WR 140 - see 
Zhekov & Skinner 2000), the exact situation is currently 
unclear with the importance of possible electron heating 
mechanisms in collisionless shocks still being debated. 

Model cw_5 lies in the middle of our grid of simula- 
tions with 77 — 0.012 and i = 45°. The latter specifies the 
orbital separation as 415/cos i — 587 AU, and the for- 
mer the mass-loss rate of the companion as a reasonable 
1.7 X lO^^Moyr"^. The computations were performed 
using grids spanning either 400 x 300 cells (cw_l — cw_6) 
or 600 X 400 cells (cw_7 — cw_8) with orbital separations 
of a minimum of 125 or 290 cells respectively. 

In Fig. 1^ we show the wind collision morphology 
for 77 = 0.02 and i = 45°. On account of its adia- 
baticity and approximately equal terminal velocities, the 
colliding winds region is stable to both thermal and 
Kelvin-Helmholtz instabilities. The intrinsic emission cal- 
culated using Raymond-Smith emissivities (with appro- 
priate abundances for the two winds) from model cw_l 
is within a factor of 2 of the measured emission (see Ta- 
ble 1^) . Since the observed luminosity possibly includes a 
significant contribution from intrinsic shocks in the wind 
of each star (cf. Owocki et al. 1988), luminosities for the 
wind-wind collision which are slightly lower than the ob- 
served value are preferred. 
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Table 2. 

1.5 X 10~ 



Parameters for the computed colliding winds models. The wind parameters for the WN8 star are M 
■Moyr-i 



and 



ASCA is 3.5 X lO-^"^ ergs 



,-1 



~ 950 kms ^ (Morris et al. 2000). The intrinsic 0.5 — 10 keV luminosity measured with 
(Skinner et al. |l999D . Since there could be additional contributions from the individual 



stars, colliding winds models with luminosities slightly lower than the observed value are preferred. Our parameter 
space study therefore incorporates this fact. Column 4 lists the inclination of our line of sight into the system (i = 0° 
is pole on, i — 90° is in the orbital plane). Column 5 and 6 list the characteristic cooling parameter for the shocked 
WN8 and OB winds. Column 7 and 8 list the intrinsic (i.e. unabsorbed) X-ray luminosity from the colliding wind 
region in the 0.01 — 10.0 keV and 0.5 — 10.0 keV bands respectively, which is dominated by the WN wind. Column 9 
lists the number of counts in the synthetic images shown in Fig. 0, assuming a 5 ksec exposure. 



Model 


n 


Mob (Mgyr-^) 


i (deg) 


XWN 


XOB 


i(o. 01-10.0) (ergs ) 


i(o.5-io.o) (ergs^^ ) 


HRC-I cts 


cw_l 


0.02 


2.8 X 10" ' 


30 


39 


2560 


2.0 X 10^'' 


4.6 X lO-"^ 


99 


cw_2 






45 


48 


3140 


1.8 X 10^^ 


3.5 X 10^^ 


76 


cw_3 






60 


67 


4430 


1.2 X 10^^ 


2.6 X 10^^ 


55 


cw_4 


0.012 


1.7 X 10"'^ 


30 


39 


4220 


0.9 X 10^^ 


2.7 X 10^^ 


56 


cw_5 






45 


48 


5170 


1.0 X 10^^ 


2.2 X 10^^ 


45 


cw_6 






60 


67 


7300 


0.7 X 10^^ 


1.5 X 10^^ 


31 


cw_7 


0.005 


7.1 X 10"® 


30 


39 


10100 


0.6 X 10^^ 


1.1 X 10^^ 


20 


cw_8 






45 


48 


12370 


0.4 X 10^^ 


0.9 X 10^^ 


17 


cw_9 






60 


67 


17470 


0.3 X 10^^ 


0.6 X 10^^ 


12 



5.^. Synthetic Images 

To help to determine if the morphology seen in Fig. || is 
the result of colliding winds emission we have calculated 
theoretical X-ray images using the hydrodynamical mod- 
els as input data. 

We would normally need to know which star is in front, 
as the stellar winds are themselves strong absorbers of X- 
rays. However, as a result of the long orbital period of 
WR 147, its orbital elements are unknown, and we do 
not know which star is in front. Nevertheless, this is not 
problematical because in WR 147 the absorption by the 
winds is actually very low (a rough estimate of the wind 
absorption along a line of sight through the WR wind to 
the region of wind collision gives A^h ~ 2 x 10^^ cm^^ cf. 



Usov 1992). The low wind column is a direct consequence 
of the wide separation of the stars and the fact that the 
sightline to the wind collision region passes through the 
outermost, least dense regions of the winds. Therefore we 
can effectively ignore the exact direction of the line of 
sight into the system, which considerably simplifies the 
calculation of synthetic images of the wind-wind collision. 

As the absorption from spherical winds is negligible, 
we assume that the total column to WR 147 (A^h = 



2.2 X 10^^ cm"^; Skinner et al. 1999) is independent of ori- 
entation. It is well known that the WR stars in Cyg 0B2 
suffer more reddening than would be normally expected. 
A correlation between stellar luminosity and reddening 
(Reddish 1967) suggests that some of the reddening is cir- 
cumstellar, possibly left over from star forma tion. For the 
extinction we use Morrison & McCammon ( 1983| ) cross- 
sections. 

As noted in the previous section, we find that the 
combined intrinsic spectrum is dominated by the shocked 
WN8 wind, so for simplicity we use WN8 emissivities for 



both winds. As Lcw,wn ^ icw,OB, this has little effect 
on the resulting emission. We also incorporate the HRC- 
I effective area from data distributed with the PIMMS 
source code, and assume a distance of 650 pc (Morris et 
al. 200C ) . The size in degrees of each cell in the hydro- 
dynamic grid varies with each model (e.g. for model cw_2 
each cell is 2.142 x 10"^ degrees square). Each HRC-I pixel 
is 3.66 X 10^^ degrees square, so for model cw_2 each HRC- 
I pixel covers the same area of sky as 292 hydrodynamic 
cells. We have therefore rebinncd the synthetic images to 
the same scale as the HRC-I pixels, and have positioned 
the WR star at coordinates (20 36 43.64, -f40 21 07.4), 
to match the centroid position of the radio contours in 
Fig. ^. The companion star has been set at a position an- 
gle of 351° from the WR-star (cf. W97 and N98). All im- 
ages are convolved with a Gaussian profile of 0.4" FWHM 
to approximate the HRC-I PSF. 

Fig. ^ shows the synthetic images from each hydrody- 
namic model assuming a 5 ksec exposure. For ease of com- 
parison, each image has the same field of view as Fig. ^ and 
the same contour levels (i.e. contours linearly spaced by 
0.366 cts). The number of counts in each image is listed in 
Table |2[ The contours of the actual observation in Fig. || 
have a small additional broadening through our use of 
CSMOOTH, which is estimated to result in a total net 
smoothing of approximately 0.42" FWHM. Therefore this 
should not overly affect the comparison. The offset men- 
tioned in Sec. || has again been applied. 

Of immediate note is the gratifying fact that the pre- 
dicted FWHM from the synthetic images (« 0.7-0.8") is 
in rough agreement with the value inferred from the data 
in Sec. |^. We also note that it appears to correlate with 
1] and i. The range in model luminosity is apparent from 
the number and extent of the contour levels in each image, 
and it is evident that model cw_l is too bright and model 
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cw_9 too faint. Because of the self-consistent generation of 
each model, the position of the peak emission lies at ap- 
proximately the same coordinates in each image. However, 
the contours appear more circular and more limb bright- 
ened for higher values of the system inclination, i (see also 
Fig. ^). For instance, models cw_3 and cw_4 have essen- 
tially the same count rate, yet the emission from model 
cw_4 is much more concentrated as witnessed by its larger 
number of contours. 

The difference in the concentration of the emission be- 
tween the models is best illustrated through profiles across 
the shock cone as shown in Fig. ||. The profile is clearly 
broader with increased rj and increased i. While this width 
can only be used in comparison with observations to de- 
termine some combination of 77 and i (they cannot be de- 
termined independently from this comparison alone) , it is 
perhaps possible for these to be uniquely determined if 
further additional comparisons between the models and 
the data are made. For instance, the normalization of the 
profile is dependent on the X-ray luminosity from the col- 
liding WN8 wind, since the shocked WR wind dominates 
the total colliding winds X-ray luminosity. With fixed val- 
ues for M vvN Bind Woown i since the shocked WN8 wind 
is largely adiabatic, Lew oc f /D. Here / is a measure of 
the shape of the wind- wind collision (and so is a function 
of 77), and the separation of the stars, -Dscp, is related to 
the seperation on the sky, Dgky, by -Dscp = -Dsky/sin i. 
Hence Lew is a (complicated) function of 77 and i. 

The width of the emission and the luminosity therefore 
yield two constraints on two unknowns, and in principle rj 
and i may be determined. We note that the observed cir- 
cularity of the emission also appears to vary, and it may 
be possible to use this as an additional constraint in the 
determination of these parameters. Due to the weakness 
of the emission from the shocked OB wind, it is, how- 
ever, clearly not possible to separately determine M ob 
and foooB from X-ray observations. 

From the contour maps in Fig. ^ we see that model 
cw_5 has 5 contour levels, in good agreement with the 
number of levels in Fig. ^ However, the number of counts 
in the synthetic image is only 45, approximately one third 
of the counts in the actual data. Although the peak bright- 
ness is approximately correct, the synthetic image is signif- 
icantly more compact than the actual data. In particular 
the synthetic emission does not extend as far south as the 
position of the WN8 star, whereas the actual X-ray emis- 
sion clearly does. For comparison, the lowest contour of 
model cw_3 extends approximately the same distance in 
RA as model cw_5, but further in DEC (+40 21 07.58 to 
+40 21 08.59). 

The intrinsic luminosity of model cw_5 is « 63% of 
the ASCA value, whereas the number of simulated HRC- 
I counts for a 5 ksec exposure is only 32% of the 
actual number of HRC-I counts. Therefore, the HRC-I 
count rate from models cw_l to cw_9 is underestimated. 
This suggests that either: i) the simulated spectrum is 



Table 3. The assumed fractional emission from each of 
the three components (wind- wind collision (WWC), WR 
wind, and OB wind) for the contour maps in Fig. |6[ 



Panel 


WWC 


WR 


OB 


left 


0.28 


0.24 


0.48 


middle 


0.28 


0.36 


0.36 


right 


0.28 


0.48 


0.24 



too hard (from comparing the collecting area of ASCA 
SIS with Chandra HRC-I as a function of energy), or 
ii) simple fits to the ASCA SIS spectum return a value 
for the characteristic absorption column that is higher 
than the actual value. These two possibilities can each 
be influenced in two ways. The spectrum could be soft- 
ened by reducing the terminal velocity of the WR star 
(since the emission from the shocked companion wind 
is insignificant its terminal velocity is irrelevant). How- 
ever, this is a fairly well defined parameter (see Morris et 
al. 200c ) . Alternatively, the spectrum could be softened if 
one or both stars contributed comparable luminosity (to 
the wind-wind collision) through emission from intrinsic 
shocks in their winds, since it is well known that such 
emission is generally softer than emission from colliding 
winds. Adressing point ii), we note that an overestimate of 
the absorbing column would be consistent with the degra- 
dation of the SIS detectors onboard ASCA (see Pittard et 
al. 200c ). One must also be careful interpreting 'character- 
istic values' when fitting single- or two-temperature mod- 
els to multi-temperature plasmas (cf. Strickland & Stevens 



1998). However, estimates from the visual extinction give 
A^H — 2.5 ± 0.4 X 10^^ cm^^, so there doesn't seem to be 
much room to maneuver in this regard. 

It is tempting to conclude from the above discussion, 
and the fact that the X-ray flux in Fig. || is more ex- 
tended in a north-south direction than predicted by any 
of the synthetic images in Fig. |[ that the most likely 
cause for this discrepency is the softening of the X-ray 
spectrum by comparable intrinsic emission from the stars. 
If true, and if the contamination is significant, it would 
seem that there is little hope that rj and i can be de- 
termined. While this would be disappointing for obvious 
reasons, we note that the S/N of the bottom data contour 
in Fig. ^ is very poor, and caution against this poten- 
tial over-interpretation. Since the absorption is so high 
towards WR 147, it is also expected that essentially all 
of the X-ray photons from the individual stars (which are 
softer than those from the wind-wind collision) will be ab- 
sorbed. On the other hand, it is indeed possible that we 
are detecting the few photons in the high energy tail of 
the stellar emission. 

We continue our analysis by performing a maximum 
likelihood fit of our models against the data. We find that 
the three models with rj — 0.02 are most favoured, and 
that the general trend is for i = 60° to also be favoured. 
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Model cw_3 appears to be preferred over the others be- 
cause its contours are simultaneously broad (as suggested 
by the East- West extension of the data in Fig. ^ and 
not too concentrated (which, for example, counts against 
model cw_l). The quality of our data does not warrant 
further detailed fitting, but there is clearly some promise 
that a future (deeper) observation could be used to per- 
form a consistency check on the best-fit parameter range 
determined by W97, N98, and CR99. 

Since we are unable to completely rule out the possi- 
bility of contamination by stellar wind photons we have 
also examined the effect of adding X-ray emission at the 
positions of the stars, where the resultant simulated image 
has again been smoothed by a suitable Gaussian. Fig. || 
shows the contours for differing ratios between the three 
components based on model cw_3, with a small adjust- 
ment in the normalization of the predicted colliding winds 
emission. Table || specifies the ratio of the individual con- 
tributions. Again the actual data and the position of the 
stars are superimposed. 

It is immediately obvious from Fig. ^ that the simu- 
lated emission is still too concentrated. However, it does 
become more extended, and in particular the addition of 
emission from the WR star spreads the contours much 
further south to give qualitative agreement with the ac- 
tual data. The addition of stellar contributions slightly 
improves the fitted likelihood in all 3 cases, although by a 
similar amount each time such that it is not clear which 
luminosity ratio is favoured. 

The biggest problem between the data and our models 
is a lack of emission in our models at declinations greater 
than approximately -1-40 21 08.6. This difference may be 
due to the incorrect alignment of the X-ray and radio data, 
the poor S/N of the lower contour levels, the emission 
from the individual stellar winds extending further from 
the stars than expected, or because our models do not 
adequately represent the actual emission from the wind 
collision region. 

4. Conclusions 

The HRC-I X-ray maps presented in this paper show evi- 
dence for the first direct (i.e. spatially resolved) detection 
of X-ray emission from the wind-wind collision zone of a 
massive early-type binary. The data from WR 147 has a 
FWHM of « 0.8", roughly double that of the PSF. Models 
of the wind collision zone are also in rough quantitative 
agreement, and in particular predict a FWHM consistent 
with that measured. To our knowledge this is the first time 
that X-ray emission has been resolved in a stellar system. 

A future exposure of longer duration is necessary to 
accurately align the X-ray and optical/radio frames, to 
potentially constrain the inclination and wind momentum 
ratio, to determine if the angular scale of the observed 
emission is consistent with models of the wind collision, 
and to determine the exact contribution of the individ- 



ual stars to the X-ray emission. Since the shocked WN8 
wind dominates the colliding winds X-ray emission, it is 



unlikely that separate values for M qb and 



be 



determined from X-ray observations of this system. 

Finally, we note the discovery that there is sometimes a 
discrepancy between estimates of the X-ray luminosity ra- 
tio of the shocked winds computed using the equations in 
Usov ( 1992D and calculations directly from hydrodynamic 
simulations. While the ease of application of the elegant 



equations in Usov (1992) have led to their widespread use, 
this issue is examined in more detail in a further paper 



(Pittard & Stevens 2002) 
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20 36 43.68 43.66 43.64 43.62 43.60 43.58 
Right Ascension (J2000) 



Fig. 2. Contours (linearly spaced by 0.366 cts) from a smoothed HRC-I image of WR 147 (thick lines) overlaid onto 
high resolution radio contours from MERLIN (thin lines). The positions of the WN8 (centroid of the thermal radio 
emission is 20 36 43.64, +40 21 07.43) and companion star (distance of 575 mas and position angle of 351° from the 
WR star) are marked with crosses. Because the X-ray image could not be accurately aligned with the optical/radio 
frame, there is some uncertainty involved in its absolute positioning. Despite this, the X-ray emission peak is certainly 
not cospatial with the WN8 star, and seems likely to be cospatial with the nonthermal radio emission from the wind 
collision zone (although it is also consistent with the position of the companion star). Here we have applied an offset of 
(—0.00965s, —0.1387") to the actual X-ray data (see Sec. ^). Statistical tests reveal the X-ray emission to be extended, 
and there is the possibility of some contaminating emission from either or both stars (see discussion in Sec. 3.4). 
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Fig. 4. Synthetic images (solid contours) of the colliding winds emission in WR147, generated from the hydrodynamic 
models listed in Table ^, and assuming a 5 ksec exposure (models are from top to bottom, left to right, cw_l to cw_9). 
The position of the WR star has been set to match the centroid of the radio contours of the southern source in Fig. || 
and the companion star has been set at a position angle of 351°. The contours are linearly spaced by 0.366 cts to match 
the levels plotted in Fig. ^. Images with more (less) than 6 contours have emission which is more (less) concentrated 
than the actual emission observed. The actual X-ray data (dashed contours) has been shifted so that its peak lines up 
with the centre of the nonthermal radio bowshock. The horizontal and vertical lines seen in some of the plots mark 
the edges of the computational domain and should be discarded as artifacts. Values for (i, rj) are given in the upper 
right corner of each image. 77 decreases left to right, and i increases from top to bottom. 
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Fig. 5. Profiles through the shock cone along the line of centres of the stars for selected synthetic images from Fig. ^ 
Each unit on the x-axis corresponds to 31 mas (i.e. 0.24 of a HRC-I pixel). The left panel displays the variation with 
i {i = 30° - solid; i = 45° - dots; i = 60° - dashes) for rj = 0.012 (models cw_4, cw_5 and cw_6 respectively). The right 
panel displays the variation with (77 — 0.02 - solid; ry ~ 0.012 - dots; r; — 0.005 - dashes) for i — 45° (models cw_2, 
cw_5 and cw_8 respectively). A deep observation which could detect these differences might constrain values of i and 
V- 
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Fig. 6. Model cw_3 with emission from the unshocked stellar winds included. In each panel the relative emission ratio 
between the three components is varied (see Table The contours are again linearly spaced by 0.366, and we have 
again superimposed the positions of the stars and shifted the actual X-ray contours. 



